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An internal singularity of a string four-dimensional black hole with second order curva- 
ture corrections is investigated. A restriction to a minimal size of a neutral black hole is 
O , obtained in the frame of the model considered. Vacuum polarization of the surrounding 

QJ , space-time caused by this minimal-size black hole is also discussed. 
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At the present time the physics of black holes contains a lot of unsolved (and even 
^ ■ non-understandcd) problems. One of them is the question on the nature of the black 

e to clarify some important aspects 
. Moreover we can also examine the 
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1 Introduction 



^^ I hole inner singularities. Studying therDLjM 

^^ ■ of the Cosmic Censorship applications Era'E 

boundaries of the applicability of the General Relativity. Another interesting and 

compktelv unsolved question is: what is the endpoint of the black hole evapora- 

[~^ I tion? B'Llfi. In order to find the more comprehensive solutions of these problems it 

0^ ' would be desirable to use the non-minimal gravity model which is the effective low 

fi I energy limit of some great unification theory. That is why during last years the 

Q^i four-dimensional dilatonic black holes attractes great attention because this type 

;' , ' £)L-h^f^^ii^r?^ni' liepi^PiSents the solution of the string theory at its low energy limit 

■ |,|^g,[|ii|iiji|izr 

It is important to note that the string theory predicts the Einstein equations 
to be modified by higher order curvature corrections in the range where the curva- 
ture of space-time has the near-Planckian values. At the present time the form of 
>^ . the higher order curvature corrections in the string effective action is not investi- 

gated completely Il3. We do not know the general structure of the expansion, and, 
hence, the direct summing up is impossible. But as we deal with the expansion, the 
most important correction is the second order curvature one which is the product 
of the Gauss-Bonnet and dilatonic terms. It increases the order of the differential 
equations till the second one and the existence of the dilatonic term makes the con- 
tribution of the Gauss-Bonnet term to be dynamic. So, the action is (for simplicity, 
only bosonian part is taken into account) 

S - ^ 

^ - 16TT 

- e-2^F^,F''^ + Afe-2^5GBJ , (1) 

Sgb = RijkiR''' — ■iRijR^-' + R ■ 
Here is the dilatonic field, F^^ = qsinOdO A dip is the Maxwell term and A is the 
string coupling constant. 
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■ThjC mo.sti.pamfuIiDvestigation of the discussed model started only a few years 
agoBEjOOE3o'E3E3li3 and all its predictions are not investigated completely yet. 
It predicts a change of the solution behavior near singularities or in the regions 
where the influence of the higher order curvature corrections becomes strong. As 
the differential equations have a very conmlicated form, the solutions were obtained 
by the perturbativeoO or numericalBliyOOE^ methods. For example, usin&theafi 
methods a new solution called "neutral Gauss-Bonnet black hole" was foundBHEil. 

The main purpose of our work is to discuss external and internal black hole 
solutions with dilatonic hair and their main properties. This means that we are 
interesting in static, spherically symmetric, asymptotically flat solutions providing 
a regular horizon. Therefore, the most convenient choice of metric (which is usually 
called as the "curvature gauge") is 

2 

ds^ = Mt^ - ^dr^ - r^{d9^ + sin^ Odip^), (2) 

where A — A(r), a = a{r). Wc use this curvature gauge (and the Einstein frame) 
for more convenient comparison with the Schwarzschild solution. 



2 Numerical results 

For searching the solution in the maximal widest range of the radial coordinate it 
was necessary to use the most "strong" method for a numerical integration of the 
systems of the differential equations with particular points. This problem was solved 
by a modernization of the methods of integration over the additional parameter Hil. 
This allowed us to investigate the internal structure and the particular points of the 
black hole using the analysis of the main determinant zeros of the linear system of 
the differential equations in the non-evident form. So, the system has the following 
matrix form 

Oil A" + ai2a' + a^cj)" = hi, (3) 

where i = 1,2,3, matrices a.y bi are: 

On = 0, 

ai2 = -m|,;cr2r-H4e~2'^A0'(cr^~3A), 

ai3 == 4e-2'/'Aa(A-CT2), 

021 = m^ia^r + 4e~2'^A</)'2Ao-, 

022 = -mlia^{A'r + 2A)-4e-'^'>'X(j)'6AA', 

023 = 4e-2'^A2AA'cr, 
asi = 4e-2'^Acr(A-cr2), 

032 = 2mlia^Ar^(t)' +Ae-^'''XA'{-3A + (T^), 

o 2 2 A 2 

033 = —Zrapicr Ar ct, 



61 = -m^;aV(0')'+4e-'^Aa(A-a2)2(,^')', 

2 r-^ 

Ae~'^'^\A(TAA'{(j}'f - 4e~2'*A0'2(A')V, 



62 = -mlia\2A' + 2Ar{^'f)-^e-^U^^ 
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Matrix 



2mli<j\4>'{A'r + 2A) - 4e-2'^A(A')V - 2e-^*q^^. 



represents the linear system of ordinary differential equations (rela- 
tively the oldest derivatives) given in a non-evident form. Hence, according to the 
existence theorem the system (0) has a single solution only in the case of its main 
determinant to be not equal to zero. In the case of zero main determinant at some 
point of the solution trajectory, the uniqueness of the solution (^) will be violated. 
Different types of zeros of the main system determinant correspond to the different 
types of the particular points of the solution and only three types of zeros present 
in the asymptotically flat solution of Eqs. (pf). So, as the main determinant has the 
following structure 



J-^niain — ^ 



AA^ + BA + C 



(4) 



where 



(-32)e-4*a2A2 



4a^6'^ml,r^ - a^ml, + Ue'^^ A' 6' \ 



B = (-32)e 



-20^4 



^A 



mpir^ + 



2e 



^20 2 



\r 



'■PI 



"^-^AV'A^ 



C 



-^^a^X^mli 



128e-6 



32e- 

+ GAe-^'^a^A'X^mlir + 128e-6'^cr^A'0'A^ 
zeros of the main determinant are 

(a) A = 0, C^O 

(b) AA^ + BA + C = 0, A^O, C^O (5) 

(c) A = 0, C = 0. 

The numerical results are shown in the Fig. 1. It presents a dependence of the 
metric functions A (a), ct (b) and the dilaton function e"^*^ (c) against the radial 
coordinate r. The solution exist in the definite range of the magnetic charge q: 
< q < rh/V2, which corresponds to the charged solution in the first order E3. The 
behavior of the solution outside the horizon has the usual form. Under the influence 
of the Gauss-Bonnet term the structure of the solution inside the horizon changes 
such that a new limiting value named "critical magnetic charge" appears. The 
structure of the inside solution is definded by this value. When qcr > q > rh/\/2 
the solution has the waiting behavior which corresponds to the fist order one. When 
< q < qcr & new singular point so-called rg appears (see Fig. 1). The type (b) of 
zero of the main system determinant is realized in this point. The solution turns 
to the other branch when reaching this point. Only two branches exist near the 
position rs- There are no any branches between r^ and the origin. Therefore a new 
singular surface with the topology S^ x R^ appears (it is an infinite "tube" in the 
t direction). Such a singularity is absent in the first order curvature gravity. This 
singularity exists in the different kinds of metric choice as we tested. 



3 The couplings r^ = rs(r,,, A) and qcr = qcr{rh, A) 

One can write an asymptotic expansions near the position r^. They are 
A = ds + diiVr - r,f + d2{Vr - rsf + ..., 



+ ox^Jr - r^ + 02{\Jr - r^) +. 
+ ^x{^r~rsf + (i>2{^r-rsf 



where r — r^ <S 1. Using these expansions and the relations between the expansion 

A2 A2(cr')2 ^ A ^AAV 



coefficients one can find the behavior of the Kretchmann scalar RijkiR^-'''^ near r^: 



RijkiR''' 
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const ^ / 1 

o 



All radial time-like and isotropic geodesies formally can be continued under the 
singularity r^ because they all have the form Const*(r — r^). Moreover, using these 
expansions one can prove that only two branches can exist near the position rg. It 
can be proved because after the manipulation with the expansions one can obtain 
only quadratic equation to the coefficient d2. 

It is possible to obtain the approximate relations r^ — r^ {r^ , A) and qcr = 

Qcr{fhi A). Substituting the Schwarzschild expansions of the metric functions (with 

the vanishing value of the dilaton charge D) into the main system determinant and 

after some manipulations one obtains {mpi — 1) 

/ \ 1/3 

r, = X'/^ 4V3r„0, (6) 



Qc 



- ^''a^^^O ■ (^) 

Fig. 2 (a) and (b) show the graphs of the values r^ and qcr against the coupling 
parameter A given by these formulas (B-H) and by the numerical integration. From 
(P) it is possible to conclude that the pure Schwarzschild solution is the limit case 
of ours with rs = 0. In the case with rather a small value of A these formulas are 
in good agreement with the results calculated by the numerical integration. While 
increasing A, the absolute error increases as a consequence of ignoring the non- 
vanishing values (for example (1 — cr), (j)'). It is necessary to point out that these 
formulas represent the dependencies r^ = const A^' '^ and Qcr = const A"'^' ^ which we 
suppose to be right because after the appropriate selection this constant by hands 
the agreement between numerical data and this formula improves. Fig. 2 shows 
also that when the influence of the GB term (or black hole mass) increases, Vg also 
increases. 



4 Minimal black hole 

4-1 Numerical investigations 

Figure 3 shows the graph of the metric function A versus the radial coordinate r at 
the different values of the event horizon r^ when q ~ 0. The curve (a) represents 
the case where r/j is rather large and is equal to 30.0 Planck unit values (P.u.v.). 
The curve (b) shows the changes in the behavior of A(r) when r^ is equal to 7.5 



P.u.v. The curve (c) represents the boundary case with r^ = rhmin where all the 
particular points merge and the internal structure disappears. The curve (d) shows 
the case where 2M <C rhmin ^-iid any horizon is absent. Here it is necessary to 
note that for minimal and near-minimal values of r^ the metric functions can be 
approximated by the following formulas 

(8) 
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- -^ 
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= Sh{rh)- 



4-2 Analytical investigations 

One can also obtain the existence of a minimal dilatonic black hole from the ana- 
lytical manipulations. This type of a system singular point corresponds to the type 
(Pa) of zeros of the main system determinant. The asymptotic expansions near the 
position rfi have the form 

A = di{r-rh) + d2{r-rhf + ..., 

a = So + si(r-r,i) + . .. , (9) 

(t) = (t)OQ + (l)iir - rh) + (l)2{r - Th)"^ + . . . , 
where {r — r^) <C 1. 

Substituting the formulas (^) into the system (^), one obtains the following 
relations between the expansion coefficients {sq, 00 = e~'^'^°° and r^ are free inde- 
pendent parameters, 

diizidj + Z2di + Z3) = 0, (10) 

where: 

Z2 = -m%irlsl, 

4 2 4 

Z3 = mpir^SQ, 
and the parameter (j)i for di 7^ is equal to: 0i = [(mp;)/(4Adi0o)] * [rhdi — Sg]. 

When di — 0, the metric function A has the double (or higher order) zero. In 
such a situation the equation for ^2 {ds, d^, . . .) is a linear algebraic one and there 
are no asymptotically flat branches. 

When di 7^ the solution of the black hole type takes place only if the dis- 
criminant of the equation (|l^) is greater or equal to zero and, hence, in this case 
T^ > 4A(/)oa/6 (q = 0, here we use the Planckian values where mpi = 1). One 
or two branches occurs and always one of them is asymptotically flat. With the 
supposition of (/)oo = (and, as we tested, in this case 1.0 < (t>Q < 2.0) the infinum 
value of the event horizon is 



C-^ = Va V4V6. (11) 

The analogous formula but in the other interpretation was studied earlier by Kanti 
et. am. In the case q > such restriction does not exist and regular horizon can 
take any meaning in the range [0 ... 00). 

Consequently, the point rhmin represents the event horizon and the singularity 
in the same point. One should remember that A is a combination of the fun- 
damental string constants. That is why formula (pT|) can be reinterpreted as a 



restriction to the minimal black hole size (mass) in the given model. This restric- 
tion appears in the second order curvature gravity and is absent in the minimal 
Einstein-Schwarzschild gravity. 

4-3 Vacuum polarization 

One of the main feathers of the black hole strnrig gravitational field i« its influence 
to the structure of the surrounding space-time BESlij'EHlLjO'Ej'ESEj. In the semi- 
classical level this effect can be described by the vacuum polarization and stress- 
energy tensor expectation values. Our purpose is to compare the < T^^, > of 
the Schwarzschild black hole with the same value of the string minimal and near- 
minimal black hole. This work is in progress now. Nevertheless, we would like 
discuss some interesting preliminary results. 

We analyzes the contribution of the massive fields to the vacuum polarization 
of the string minimal black hole. In this case the contribution will be considerable 
enough because its mass has the order of the Planck mass mpi = yJficjG. Work- 
ing almost near the boundaries of applicability, we can study the expansions for 
Hartle-Hawkiag vacuum average of the stress-energy tensor. Using Frolov-Zel'nikov 
expansions Erlla 

<T(^)>--^^^ (12) 

which in one-loop approximation reads 

W^il = y-^^/d^x|5r/^L(^)+0(a, (13) 

(4:TTm)^ 8 7'. J 

-^ — O-lRafi-ia-eR'^ ''"^'^ + a2R'pRa-/creR ^'^'^ + O.sR'l'JjRt^^R^sRatS 

+ asRa/3--fR + agR^RpR^ + aioRa/sRjaR ■ 

In our coordinates (A, cr, see expansions (^)) the stress-energy tensor average 
values (|lj) have the following form (because of the great size of the formulas we 
show only the < Tq > component in the case when the spin of a particle is equal 
to 1/2) 

< TqO > = T1/T2, (14) 

- ^^^r'^^rlsh - 17962560r^'^Sft + 108404640r^3 , 

31 

- 179223493r^V^s,, + 89133599r^V^Sft - ^ 

oX 

+ ^-^r^'r.sl - ^-^r^'rlsl - 583036720r«S^, 

+ 2813006280rSV,,s^ - 4012428650r'5V^s^ + 1783041602r^V^s^ 
32928 fin . 567518 ,„ , 528220 ,0 n , 
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6468 68 2 



- 2729994080r^®s}^ + 11437564864r5VftS^ - 14744667622r 
+ 6035301470r^^'V^^ - ^-^r''4 + ^^-^^,4 

- ^^^r^"rlsi - 2818841780r5"4 + 10549836302r^'> r^s^ 

74039 44 5 



,56„2 3 



12479930891r''*r^4 + 4749439473r''V^s;^ + ^^ 

405230 43 5 359660 42 2 5 r.noroo.o 42 5 

^^^ ''^/.^ + -^^r^ rlsl - 540258212r42s5 

1847352418r''V,,s^ - 2027148281r4"r^s^ + 720283059r^^r|s^ 
46109 36 6 98490 35 6 , 115640 34 .^ 

3278996r34s6 - 4899564r33rftS^ + 1843020r32r^s| 
18963 28 7 , 86534 2, , 87612 2 2,7 



i>i di di 

1820140r26s^ + 3130052r25r,,s^ - 1456108r2V^s^ 

Zg^,20,8 _ ^,i9,,4 ^ 47i24ri«4 
24438riV,.« + ^r^^.^ + |r"r,.)^ 
21148ri°4 - 2790r9r^s^ - —r\s]^ 



T2 = r88 - llr8"s^ + 55r^2s2 _ ^g5^64^3 

+ 330r''^S^ - 462r'^®s^ + 462r''°s^ - 330r3 
+ 165r2V, -55ri6s? + llr«5i°-s". 



< T.^ >^4) = 



AsP^^ + 3i2Bs + As)qi^ 



Here one can see that when Sh = 0, then < Tq > takes its Schwarzschild value 

'" ^ 100807r2m2r8 
where 

pH -15 + 16^ p;^-3 + 42M p0^p^^Q_n2M 

ql = -10 +11^ q;^4-3^ ql = q^^~l2+lA^ 
Comparing our stress-energy tensor with the Schwarzschild one, we see that 
their numerical values coincide asymptotically at the infinity, but strongly differs 
in the neighborhood of the event horizon (for the Planckian masses the difference 
is about 15 orders). Therefore taking into account the terms of the string theory 
shows its influence to be stronger near the event horizon where the gravitational 
field is not a week one. So, we arrive to a conclusion that the application of the 
General Relativity gives a good results at the large distances from this minimal 
black hole, but in the neighborhood of this object one must use quantum gravity 
models. 



5 Discussion and conclusions 

From all the discussed numerical and analytical arguments it is necessary to con- 
clude that string theory gives the new knowledge on the black hole space-times. It 
provides new types of the internal singularities which are absent in the same models 
of the General Relativity. 

The most interesting consequence of the second order curvature corrections 
is the existence of a minimal black hole. It exists-,only in the neutral (quasi- 
Schwarzschild) case. Using Green- Witten-Schwartz E3 formulas (Chapter 13), we 
can calculate the value of the string coupling constant A, therefore, we can find the 
numerical value of this minimal black hole mass. It has the value about 0.4 mpi 
(and, in that case, st = 1/40). Speculating about this phenomenon (if this object 
is stable) we can say that in General Relativity the most realistic models describing 
theJaJi 
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:k holes in our Universe are Schwarzschild one and Kerr one. Sometimes 
it was established that spinning black hole looses its angular momentum 
during the Hawking radiation. Therefore, our minimal (quasi-Schwarzschild) black 
hole being the endpoint of the evaporation can represent the relic remnant of the 
black holes formed during the initial stages of our Universe formation. This is the 
very interesting problem and it requires the additional investigations. 
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Figure 1: 



Dependence of the metric functions A (a), a (b) and the dilaton function e ^"^ 
against the radial coordinate r when < g < qcr 



(c) 



Figure 2: 



The dependencies o 
the formulas M 



.the values rs and Qcr against the coupling parameter A given by 
(lines) and by the numerical integration (squares) 



Figure 3: The dependence of the metric function A versus the radial coordinate r at the different 
values of the event horizon r;, when q = 0. The curve (a) represents the case where r/j is rather 
large and is equal to 30.0 Planck unit values (P.u.v.). The curve (b) shows the changes in the 
behavior of A(r) when r^ is equal to 7.5 P.u.v. The curve (c) represents the boundary case with 
^h = ^hmin where all the particular points merge and the internal structure disappears. The curve 
(d) shows the case where 2M <g r^^^^„ and any horizon is absent 



Figure 4: The dependencies of | < Tq > (r)| at string case with the second order curvature 
corrections (a) and at Schwarzschild (b) case 
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